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ABSTRACT 

Recent astronomical data strongly suggest that a significant part of the dark matter, 
composing the Local Group and Virgo Supercluster, is not incorporated into the galaxy 
haloes and forms diffuse components of these galaxy clusters. Apparently, a portion 
of the particles from these components may penetrate into the Milky Way and make 
an extragalactic contribution to the total dark matter containment of our Galaxy. 

We find that the particles of the diffuse component of the Local Group are apt to 
contribute ~ 12% to the total dark matter density near the Earth. The particles of the 
extragalactic dark matter stand out because of their high speed (~ 600 km/s), i.e. they 
are much faster than the galactic dark matter. In addition, their speed distribution is 
very narrow (~ 20 km/s). The particles have isotropic velocity distribution (perhaps, 
in contrast to the galactic dark matter) . The extragalactic dark matter should give a 
significant contribution to the direct detection signal. If the detector is sensitive only 
to the fast particles (v < 450 km/s), the signal may even dominate. 

The density of other possible types of the extragalactic dark matter (for instance, 
of the diffuse component of the Virgo Supercluster) should be relatively small and 
comparable with the average dark matter density of the Universe. However, these 
particles can generate anomaly high energy collisions in direct dark matter detectors. 

Key words: dark matter, elementary particles, galaxies: haloes, solar neighbourhood. 



It is widely believed that all the dark matter particles (here- 
after DMPs), which a terrestrial observer can detect, belong 
to the Milky Way Galaxy. The main aim of this letter is to 
dispute this assertion and to show that a remarkable frac- 
tion of dark matter particles detected on the Earth does 
not probably belong to our Galaxy. Although their density 
is relatively small, as compared with the total dark matter 
density, their impact into the direct detection signal may 
even dominate because of high speeds of the particles. 

According to the modern cosmological notion, the 
haloes of giant galaxies, like Milky Way or Andromeda, 
are regions of local dark matter overdensity, rather than 
isolated islands. Indeed, the Local Group, along with the 
haloes of large and dwarf galaxies, contains a significant 
fraction of dark matter that is not bound in the galax- 
ies and presumably forms a larg e envelope of the Local 
Group (|Binnev fc Tremainel [2008). A significant part of 
the dark matter of the Virgo Supercluster is also not lo- 
calized in haloes and probably distributed more or less 
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homogeneously over all the vo lume of the Supercluster 
(Mak arov fc Karachentse"vll201ll ). Some part of this diffuse 
dark matter (preeminently from the Local Group envelope) 
penetrates into the central region of our Galaxy and gives a 
contribution to the direct detection signal, which can even 
dominate under certain conditions: as we will see, the den- 
sity fraction of the extragalactic dark matter is relatively 
small (~ 12%), however its particles should have extremely 
high speeds, close to the escape velocity (~ 600 km/s) or 
even higher. It sets off the extragalactic particles from the 
halo DMPs with much lower average speed. The direct de- 
tection signal produced by this component should also have 
some other characteristic features that will be discussed be- 
low. 



2 THE DARK MATTER ENVELOPE OF THE 
LOCAL GROUP 

Unfortunately, the total mass, distribution and dynamical 
properties of the extragalactic dark matter environment are 
now poorly known. Therefore, we have to do with rough 
estimates of its content near the Solar System. 
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The Local Group consists of two very massive galax- 
ies (Milky Way and Andromeda galaxy M31), less mas- 
sive Triangulum galaxy M33, and a host of dwarf galax- 
ies. It seems reasonable to say that the Local Group con- 
tains a massive diffus ed dark matter component as well 
<|Kahn fc Woltiej|l959h . Unfortunately, some parameters of 
the system have not been measured with the adequate ac- 
curacy. We accept the following values in this letter: the 
radius of the Solar System orbit Iq = 8 kpc, the Milky 
Way mass Mmw = 10 12 Mq, the Andromeda galaxy mass 
M31 = 1.6xl0 12 M( 7) , the distance between them d — 750 kpc 
l|Cox fc Loebl [20081 ) . The tangential components of the ve- 
locities of even some massive members of the Local Group 
are also not quite explored, and we know almost nothing 
about the distribution and dynamical parameters of the dif- 
fused component: so the investigation of its motion in the 
complex gravitational field of several bodies is quite a diffi- 
cult and underdefined task. However, our aim is much sim- 
pler: we would like to model just the process of the envelope 
dark matter penetration towards the Solar System. Taking 
into account all the above-mentioned uncertainties, we will 
try to construct a simple toy model, that does not claim to 
describe all properties of the Local Group, allowing us, how- 
ever, to estimate the density and the velocity distribution of 
the extragalactic dark matter near the Solar System. 

Let us consider the following model: the system is sta- 
tionary and spherically-symmetric. Specific angular momen- 
tum fi = [v X r\ and maximum radius ro the particle moves 
from the centre remain constant for each particle in such 
a system, and gravitational potential 4>(r) depends only on 
r. ro of the particles belonging to the envelope lie in some 
interval [r in ,r ou t\. We accept r in = 300 kpc, which approx- 
imately corresponds to the size of the Milky Way Roche 
lobe in system Milky Way - M3 1, r ou t = 600 kpc, in ac- 
cordance with (|Cox fc Lo cb 2008). The particles have some 
distribution /(ro) over ro inside [r; n ,r ollt ]; we assume that 
their specific angular momentum fj, = [v X f] has Gaussian 
distribution. So the overall distribution (i.e. the mass dm of 
particles in some interval drod/x) is: 



dm = /(ro)-=2 exp 



2 ± 



dfidro, 



r G [ri n , rout] (1) 



where a is, generally speaking, a function of ro. We accept 
the envelope mass M env — J r " ut f(ro)dro = 1O 12 M0. This 
value is noticeably smaller than the total mass of the diffused 
compone nt, which is estim ated as ~ Mmw + M31 = 2.6 x 
10 12 M© (|Cox fc Loebl200Sh . We allow for, however, that the 
main part of this substance surrounds and accretes on the 
Andromeda galaxy, and take only Mmw / (Mmw + M31) of 
the total mass, so we assume that the mass should be divided 
proportionally to the Roche lobe areas of the components. 

At first glance it would seem that the above-stated 
model is completely unusable to describe the dark matter 
motion in the Local Group: the gravitational field of the 
system by no means can be considered as central on the 
scale ~ 600 kpc, because of huge perturbations from other 
group members, especially from the Andromeda galaxy. 
However, we are only interested in the envelope dark matter 
penetration towards the Solar System, and this process is 
totally defined by the angular momenta of the particles: only 
the particles with very small momenta can reach the Earth. 
The motion of the particles near the Solar System is almost 



unaffected by M31 (Iq <C d) and may well be described by 
the above-mentioned model. When a particle moves from 
the envelope to the Earth, its angular momentum is, of 
course, strongly influenced by the tidal perturbations. How- 
ever, we almost do not know the momentum distribution of 
the particles in the envelope. Therefore equation ((TJ) may 
be thought of as describing the resultant distribution of the 
falling particles with regard to the perturbations from the 
other members of the Local Group. Moreover, as we will 
see, the shape of angular momentum distribution is not very 
important: we actually use only the value at [i = 0. As for 
perturbations of the particle energy, they are of the order of 
GMzi/d, i.e. always small. This is a result of the fact that the 
main part of the particle acceleration takes place deep in the 
Milky Way, where the gravitational field is much stronger. 

Now we should find the particle distribution inside ra- 
diu s Tin. A very s imilar task has been studied extensively 
m (|Baushevll2012ar i. We cite here only the results adaptable 
to our work, skipping the complete derivation. The exact 
distribution inside n n depends on r and is equal to: 



P = 



/(r )roAtexp (-/i 2 /« 2 ) d^dr 



27rra 2 (r )T(r ,^)\/rg -r 2 yfl 



(2) 



Here T(ro,fi) is the half-period of a particle with maximal 
radius ro and specific angular momentum fj,, i.e. the time it 
takes for the particle to fall from its maximal radius to the 
minimal one, and fi ma x is the maximum angular momentum 
of a particle wherewith it can reach radius r 



2 

Pmax 



2(0(ro) 



""•l<(3-3 
'0 



(3) 



Our concern is only with the particle distribution at r = Iq . 
Since Iq <C r» n < ro, we can simplify the above equations 



P = 



: /(ro)/xexp (— /x 2 /a 2 ) d^dr 



27r^0Q 2 (ro)T(ro)v74^r _— M 2 



(4) 



(•>) 



Here we took into account that for /1 £ [0, ^ m ax(lQ)] period 
T(ro,fi) is almost independ ent on p (T(ro,/x) ~ T(ro,0) = 
T(r ), see (|Baushevll2012ah for details). 

Equation (J4j can be significantly simplified, if we take 
into account that a(ro) is physically constrained. On the one 
hand, a(r ut) hardly can be higher, than 



a(r otl i 



1G(Mmw + M env ) 
r 'out 



(6) 



since in the opposite case a significant fraction of the 
particles would have the speed above the escape velocity 
at this radius. On the other hand, numerical simulations 
jStadel et al.| [2009) show that the root-mean-square angu- 
lar momentum of the particles should be quite high and close 
to upper limit © . Though there is rather strong evidence 
that a of the particles of our Galaxy is much lower (|Baushevl 
l201ll) . we will use value ([6} in our calculations, since the den- 
sity of the extragalactic dark matter grows with decreasing 
of a, and we take the maximum possible value in order to 
obtain a conservative estimate. The dependence of a on ro 



© 2012 RAS, MNRAS 000, \TM 



Extragalactic dark matter and direct detection 3 



is not well known, and we will presume it to be power-law 
a(r ) = a(r out ) (j-O (7) 

\ Tout J 

However, we only need much softer condition a(ro) ^ 
fJmaxQe) ^° simplify ©• It means that the tangential ve- 
locity dispersion in the envelope is supposed to be higher 
than « esc ^- ~ 16 km/s. Such an assumption seems quite 
natural. If a(ro) ^ Pmax(lo), we can simplify Q as 



f(ro) fJ.dfid.ro 



o 2nl e ,a 2 (r )T(r ) \/ pliax 



(8) 



Now we should ascertain the velocity distribution of the 
particles. Let us denote the tangential and radial compo- 
nents, and the total velocity of a particle at r — Iq by u T , 
u r , and u respectively. We can use u T and u instead of ro, 
fj. Indeed, fi = u T l Q , u = ^/2(0(r o ) - <j>(l®)), fJ max = ul®. 
Let us consider the particles with the same u (i.e., with the 
same ro) and find their angular distribution. An element dQ 
of the solid angle in the phase space is equal to 



dtt = 



U 1"KU T du T 



fidfi 



u r 4tt« 2 2/q ^fl 2 max - p 2 v/2(0(r o ) - 0(/ )) 

(9) 

Substituting this equation to ((8}, we obtain 



P = 



/(ro) x /2(0(ro)-0fe))dr o 
7TQ 2 (ro)T(ro) 



dtt 



(10) 



As we can see, the particle distribution depends only on ro 
(i.e., on velocity magnitude u), and is independent on the 
direction. So the distribution is isotropic, u and ro are bound 
by a one-to-one relation 

, (^f)dr 

u = Vmro) - 0(i©)) du = ) ° L == (11) 

V2(<^(ro) - W®)) 

We can substitute this equation to (|10p and take into ac- 
count that J dQ. = 47r. 



P = 



8/(r o )(^(r o )-0fe)) 
a 2 (r )T(ro)(#(ro)/dro 



■du (12) 



% /2(^.(r I „)-0(I )) 



Hence the velocity distribution of the dark matter from the 
envelope is isotropic near the Solar System; equations (|11[) 
and (|12[1 totally define its density and momentum distribu- 
tion. 

In order to complete the solution, we should define func- 
tions 4>{ro), T(ro), and /(ro). We will assume /(ro) to be a 
power-law function with some index j. Since J /(ro)dro = 

j^tenv i 



f(ro) = 



(j + l)M env 



T out Ti n \ Tout Ti; 

We accept the mass distribution at r > r; n to be 
M(r) ~ M w + M en „ ' 



(13) 



(14) 



v Tout Ti; 

This equation is not quite true, because in fact the particles 
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Figure 1. The normalized velocity distribution of dark matter 
particles near the Earth in the frame of reference that does not 
rotate around the Galaxy centre. The distribution of the galactic 
DMPs is supposed to be Maxwell 1231 1. The extragalactic compo- 
nent gives a narrow high peak near 600 km/s. 



contribute as well to the mass inside r; n . However, our cal- 
culation is estimative, and we may disregard rather a small 
deviations from (|14|) . With the same accuracy T(ro) may be 
thought of as a power-law function. One can readily see that 
T(ro) oc r 1_ 5 if M(r) oc r 3+1 . Though the mass distribution 
is not quite power-law in our case, we presume 



T(r ) = T(r in ) I H 



(15) 



and T{ri n ) — 10 17 s, which is slightly more, than the time 
necessary for a body to fall on point mass 1O 12 M0 from 
300 kpc with no initial velocity. Function 

is totally defined by (111ft . 



dr 



GM(r ) 
3 

1) 



#(ro) = G 
dr rl 



Mmw + M env 



ro 



Tout 



(16) 



As for function y/2(cj>(ro) — 4>(Iq)), it remains almost con- 
stant for ro € [r in ; r out ] owing to the smallness of \4>(r out ) — 
4>{rin)\ as compared with \(j>(ri„) — 0(r; )j . Therefore we can 

approximate ^/2(tj>(r ) - <£('©)) ^ \/2(<j)(rin) - 0(Z©)) = 
V. Hereafter we accept V = 600 km/s. The velocities of 
all the particles lie in a very narrow interval AV 

AV = ^V 2 + 2{<t>{r ut) - W^nf) ~ V (17) 

Now we should substitute equations 0, (|15|) . (|13|l for a(ro), 
T(ro), /(ro) respectively and (O for a(r ou t) into equation 
(|12|l . It is convenient to introduce k = r out /r in — 2. After 
some trivial calculations we obtain 



p = 9 



+ k 2j (ki j - 2i - 1) 



V 



(|j-2i) (fc-l)^ 1 GT(r ln )r ut 



(18) 



The shape of velocity distribution function n in the phase 
space is given by 



n(u) oc 



du 3 

4W'° 



|j-2i-i ( dro 
du 



(19) 



where we should substitute ro by u G [V, V + AV] with 
the help of equations {13]| and (fl6)l . Equations (Tl8)) . (fl9)l 
completely determine the solution of the task. 
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3 RESULTS AND DISCUSSION 

Not much is known about the dark matter distribution 
in the envelope. On the other hand, as we can see from 
equation (|18|l . the result is not strongly dependent on the 
choice of i and j (as a consequence of the relative small- 
ness of ratio r out /ri n ). It seems reasonable to choose i and 
j by analogy with the well-known isothermal halo solution 
(dM/dr = const and the Maxwell DMP velocity distribu- 
tion with a temperature, constant over the halo), which 
corresponds to i = 1, j = 0. Substituting these values in 
combination with T(ri n ), a(ri n ), and V to (|18[1 . we obtain 
the density of the extragalactic dark matter near the Earth 
p — 3.7 x 10 -2 GeV/cm 3 . The speed distribution is no- 
tably narrow: the absolute values of all particles fall within 
AV ~ 16 km/s. As we have already mentioned, this fea- 
ture is a consequence of the smallness of \4>{r out ) — 4>{ri n )\ 
as compared with —4>(Iq). Therefore, two properties of the 
velocity distribution are model-independent: the speeds of 
extragalactic DMPs from the envelope lie in a narrow range, 
and their angular distribution is isotropic. 

The density of the extragalactic dark matter turns out 
to be fairly high: 3.7 x 1(T 2 GeV/cm 3 is more than 12% of 
the total dark matter density near the Earth ~ 0.3 GeV/cm 3 
IIGorbunov fc Rubakov|[201ll ). This brings up a question: 
How reliable is the estimation? Above we have already dis- 
cussed the approximation of the system by a spherically 
symmetric model and found it acceptable. The premise that 
/(ro) terminates abruptly at ri n and r ou t is also unphysical. 
Undoubtedly, our result is assessed; however, it cannot be 
called optimistic. Indeed, as the dimensional method shows, 
for any envelope model the density of the extragalactic dark 
matter is, with an accuracy of a numerical factor, equal to 

p0C l^W (20) 

where (a) and (T) are the average values of the respective 
quantities. Our calculations confirm this dependence: it can 
be easily obtained from \V2). v esc is almost independent on 
the model choice. (T) is essentially defined by the size of the 
Milky Way Roche lobe, and thus is also more or less model- 
independent. The main source of the uncertaint y is envelope 
mass M env . We proceed from the assumption of lcox fc Loebl 
(2008) that M env is approximately equal to the masses of the 
galaxies of the Local Group. We used the highest possible 
value (JS| for a: it a was higher, a significant part of the 
envelope would rapidly evaporate. Since p tx a -2 , this choice 
is conservative. 

Thus there are two possible situations, when (|18|1 signif- 
icantly overestimates the density of the extragalactic dark 
matter. It may be so, if the envelope mass is in fact much 
lower than the masses of the Local Group member galaxies. 
The strong overestimation may also appear, if the angular 
momentum distribution of the envelope DMPs differs greatly 
from the Gaussian {TJ, i.e., almost all the particles have cir- 
cular orbits. Such a supposition seems highly improbable. 
First of all, it is in sharp contrast to N-body simulation 
results l|Stadel et al.l 120091 ) . There is also a good indirect 
counterargument: the largest Local Group member M31 has 
quite low angular momentum and, consequently, very oblong 
orbit dKahn fc Woltierlll959l '). It is plausible that Milky Way 
and M31 will finally experience a central collision. Thus the 
presence in the diffuse component of the Local Group of a 
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Figure 2. The ratio of the direct detection signal produced by 
the mixture ~ 12.3% of extragalactic component and 87.7% of 
the galactic DMPs to the signal produced by pure galactic dark 
matter. The fraction and velocity distribution of the extragalactic 
component were calculated in accordance with (1 1 S D - We consid- 
ered two models of the velocity distribution of the galactic DMPs: 
Maxwell {23) (solid line) and anisotropic | |24[ ) (dashed line). The 
extragalactic dark matter almost does not affect the signal, if 
Vmin < 300 km/s, but totally dominates above 450 — 500 km/s. 



bulk of dark matter particles that have very oblong orbits 
and can reach the Earth seems quite possible. 

The 12% extragalactic component of the total dark mat- 
ter density can be especially important for the direct dark 
matter search. The direct search is based on the detection 
of the collisions of dark matter particles with nuclei of the 
target. The signal is sensitive to the velocity distrib ution: 
roughly speaking Eela nger. Nezri. fc Pukhovl (|2009l ). it is 
proportional to 



^ld 3 v 



(21) 



Here fi(v) is the distribution in the Earth's frame of ref- 
erence: it should be obtained from (O, (f2"3")l . or (f2H) 
by a Galilean transformation. v m i n is the minimal DMP 
speed, to which the detector is sensitive (see details in 
iBelanger. Nezri. fc Pukhovl (|2009T l'). 



E A {m x + m A ) 
2 mAirii 



(22) 



where m x and tua are the DMP and the detector nucleus 
mass respectively, E A is the detector activation energy, de- 
pending on its construction. In order to estimate the influ- 
ence of the extragalactic dark matter to the direct detection 
signal, we should define a model for the velocity distribution 
of the galactic component. The Maxwell distribution is now 
routinely used, mainly because of its simplicity: 



n{u) 



(0Fv ) 3 



exp 



(23) 



vq is the orbital speed of the Solar System. There are strong 
reasons to suppose, however, that the distribution of the 
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galactic DMPs is strongly anisotropic and looks like 



exp 



n(u) — 



2ol 



1 

max 



(24) 



where u r £ [~v max ;v max ], u max — 560 km/s, ao = 80 km/s, 
u r and u T are radial and t angential comp onents of the par- 
ticle velocity, respectively |Baushevll201ll '). 

We calculated the signal produced by mixture of 
3.7 x 10 -2 GeV/cm 3 of the extragalactic dark matter and 
0.263 GeV/cm 3 of the galactic one and divided it by the 
signal produced by the pure galactic dark matter with the 
same total density (0.3 GeV/cm 3 ). We used both models of 



the galactic DMP distribution: Maxwell ()23[) and anisotropic 
(|24[) . The two ratios are represented in Fig. [2] by the solid 
line (Maxwell velocity distribution of the galactic dark mat- 
ter particles) and by the dashed line (anisotropic velocity 
distribution). One can see that the signal is scarcely affected 
by the presence of the extragalactic component, if v max < 
300 km/s. However, the situation drastically changes for 
higher v max ; if v max is larger, than 450—500 km/s, the extra- 
galactic signal dominates. This is not particularly surprising: 
all the extragalactic particles are faster than 600 km/s, while 
the number of the galactic particles rapidly drops above 
~ 450 km/s. Hence the impact of the extragalactic dark mat- 
ter to the direct detection signal can be very important, es- 
pecially if the DMP mass m x is small. Indeed, if m x is small, 
v max is high ()22|) . i.e., we can detect only th e fastest DMPs. 
For i nstance, DAMA collaboration reports (Be rnabei et al.l 
120111 ) about the detection of a signal produced by ~ 10 GeV 
weakly interacting massive particles. We shall not discuss 
here the question of the natu re of the signal (o ther detec- 
tors do not confirm the result l|Aprile et alJl201lh ). It should 
be recorded, however, that v max > 450 km/s for the major- 
ity of detectors, if the DMP is so light, i.e., the extragalactic 
component signal should totally dominate. 

Notice that we should, strictly speaking, have cut dis- 
tributions (|23|) and (|24l) at u — v esc . However, the fraction 
of the particles with u > v e3C is negligible for both the dis- 
tributions, and the cutting would hardly affect the result; 
the impact of the extragalactic component would be even 
slightly higher. 

In conclusion, let us briefly consider the extra- 
galactic dark matter that does not belong to the Lo- 
cal Group. As recent astrono mical observations imply 
ijMakarov fc Karachentsevl l201lT l. the dark matter of the 
Virgo Supercluster, in addition to galaxies and their groups, 
forms a large diffuse component. We do not know its ve- 
locity and space distributions, but it seems reasonable to 
assume that the dark matter is distributed more or less uni- 
formly, and the velocity dispersion of the DMPs is compara- 
ble with that of the observable members of the Supercluster 
(voo ~ 500 km/s). The measurements estimate the average 
density of the diffuse component as p ~ 10 -6 GeV/cm 3 . The 
gravitational field of the Local Group should increase this 
quantity near the Earth. However, we can roughly estimate 
the enhancement as 1 + v^^/v^, where v esc ~ 65 km/s is 
the es cape velocity from the Solar System orbit |Baushevl 
l2012bl ). Thus the density of the Supercluster dark matter 
is approximately 3 times higher near the Solar System, but 
yet hardly exceeds 10 -5 GeV/cm 3 . This value is so low, 
that it may scarcely be of interest for modern experiments. 



On the other hand, the Supercluster DMPs are particularly 
energetic (v > 1000 km/s) and hence may give a very char- 
acteristic signal. 
To summarize: 

1) The particles of the diffuse component of the Local 
Group are apt to contribute > 10% to the total dark matter 
density near the Earth. 

2) The particle speeds are ~ 600 km/s, i.e. they are 
much faster than the galactic DMPs. The particles have 
isotropic velocity distribution (perhaps, in contrast to the 
galactic dark matter); their speed distribution is very nar- 
row (AV ~ 20 km/s). 

3) The extragalactic dark matter should give a signif- 
icant contribution to the direct detection signal. If the de- 
tector is sensitive only to the fast particles (v > 450 km/s), 
the signal may even dominate. 

4) The density of other types of the extragalactic dark 
matter (for instance, of the DMPs forming the diffuse com- 
ponent of the Virgo Supercluster) should be relatively small 
and comparable with the average dark matter density of 
the Universe. However, these particles can generate anomaly 
high-energy collisions in direct dark matter detectors. 
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for the contents of this publication. 
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